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To quantify the reduction in workload during intra-aortie balloon pump (IABP) therapy, 
indirect parameters are used, such as the mean arterial pressure during diastole, product of 
heart rate and peak systolic pressure, and pressure-volume area. Therefore, we investigated 
the cardiac energy consumption during IABP therapy using a cardiac electromechanics 
model. We incorporated an IABP function into a previously developed electromechanical 
model of the ventricle with a lumped model of the circulatory system and investigated the 
cardiac energy consumption at different IABP inflation volumes. When the IABP was used 
at inflation level 5, the cardiac output and stroke volume increased 11%, the ejection 
fraction increased 21%, the stroke work decreased 1%, the mean arterial pressure 
increased 10%, and the ATP consumption decreased 1 2%. These results show that 
although the ATP consumption is decreased significantly, stroke work is decreased only 
slightly, which indicates that the IABP helps the failed ventricle to pump blood efficiently. 

Key Words: Intra-Aortic Balloon Pump; Cardiac Electromechanics Model; ATP Consumption; 
Stroke Work 



INTRODUCTION 

The intra-aortic balloon pump (IABP) is used to increase myo- 
cardial oxygen perfusion, while at the same time increasing car- 
diac output and decreasing the workload of the ventricle. This 
is realized via counter-pulsation of the IABP. Several parameters 
affect IABP treatment, such as the volume of the balloon, its lo- 
cation in the aorta, its rate of inflation and deflation, and the 
phase of IABP pumping. Several researchers have examined 
the effects of various parameters on IABP function in terms of 
increasing myocardial oxygen perfusion and decreasing the 
workload of the ventricle. When the volume of blood used to 
inflate and deflate the balloon increases, the IABP functions 
were improved (1). Inflation of the IABP at end-systole and de- 
flation of the IABP at end- diastole were the optimal pumping 
phase in order to decrease the ventricle workload and increase 
aortic pressure during diastole, thereby increasing myocardial 



oxygen perfusion (2). 

To quantify the reduction in workload, indirect parameters 
are used, such as the mean systolic pressure (3), product of heart 
rate and peak systolic pressure (3), and pressure-volume area (4). 
No study has quantified the reduction in workload with IABP 
directly because experimental methods for documenting and 
evaluating myocardial energy consumption throughout the 
ventricular volume are hampered by low spatiotemporal reso- 
lution. Computational modeling is an alternative approach that 
overcomes this limitation. However, the previous computation- 
al models of the cardiovascular system and IABP used a simple 
lumped model of the cardiovascular system, which cannot cal- 
culate the contractile cardiac energy consumption directly. 

Previously, we developed a computational model of left ven- 
tricular assist device support using a three-dimensional elec- 
tromechanical model of failing canine ventricles together with 
a lumped model of the circulatory system (5, 6). This model en- 
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abled us to quantify the changes in the local contractile energy 
consumption of the myocardium. In this study, we incorporat- 
ed IABP function into the electromechanical model of a ventri- 
cle with a lumped model of the circulatory system. Then, we 
used the model to investigate the changes in the contractile en- 
ergy consumption of the myocardium with an IABP using dif- 
ferent inflation volumes. 

MATERIALS AND METHODS 

Model description 

The IABP-implanted cardiovascular system model was imple- 
mented by combining a three-dimensional (3D) finite element 
(FE) electromechanical model of the failed heart, a lumped mod- 
el of the circulatory system, and a model of IABP function. Fig. 1 
shows a schematic of the integrated model (A) and myofilament 
dynamics, which were adapted from Rice et al. (7). The models 
are described in detail below. 

3D cardiac electromechanical model of the failing heart 

We used an anatomically accurate 3D-FE electromechanical 
model of failing canine ventricles, which incorporated a bio- 
physical representation of cardiac myofilament dynamics (8). 
The electrical model based on the membrane dynamic model 
of Fox et al. (9) used magnetic resonance (MR) image-based 
canine ventricle geometry with approximately 1,400,000 nodes 
and 1,700,000 elements, which were developed using linear in- 
terpolation. The mechanical model based on the myofilament 



dynamic model of Rice et al. (7) used the same canine ventricle 
geometry with 356 nodes and 172 elements, which were devel- 
oped using Hermite polynomials. Both meshes also have realis- 
tic fiber and laminar sheet structural information obtained from 
the diffusion tensor (DT) MR image. The passive mechanical 
properties of the myocardium were assumed to be transversely 
isotropic, hyperelastic, and nearly incompressible, defined by 
an exponential strain energy function. The local electrical acti- 
vation times in the electrical mesh were mapped to the mechan- 
ical mesh. Fig. 2 shows the electrical activation time mapped to 
the mechanical mesh. Here, the local electrical activation time 
is defined as the instant at which transmembrane voltage ex- 
ceeds 0 mV 

As an electrical wave propagates through the heart, the de- 
polarization of each myocyte initiates the release of calcium (Ca) 
from intracellular stores. This is followed by the binding of Ca to 
troponin C and cross-bridge cycling. The cross-bridge cycling 
forms the basis for contractile protein movement and the de- 
velopment of active tension in the cell, resulting in deformation 
of the ventricles. These physiological contractile mechanisms 
are considered in the myofilament model (7). 

To consider remodeling of the passive electromechanical 
properties associated with heart failure (HF), several changes 
were made to the electromechanical model (Table 1). First, we 
used the dilated geometry of a canine ventricle with heart fail- 
ure (HF) (8). Then, the electrical conductivity was reduced by 
30%, allowing for a total electrical activation time of 150 ms, ac- 
cording to the experimental results of Helm et al. (10). Third, to 
account for the increased stiffness of the failing myocardium, 
the passive scaling constant in the strain-energy function was 
increased fivefold (11). Finally, to introduce systolic dysfunction, 




C S A,iABp(t) = CsA-ScaleFactoriABp(t) 



Fig. 1. Schematic diagram of the finite-element ventricular electromechanical model 
coupled with the circulatory model (A). Prv, RV pressure; Vrv, RV volume; Plv, LV pres- 
sure; Vlv, LV volume; Rpa, pulmonary artery resistance; Cpa, pulmonary artery compli- 
ance; Rpv, pulmonary vein resistance; Cpv, pulmonary vein compliance; Rmi, mitral 
valve resistance; Cm, left atrium compliance; Rao, aortic valve resistance; Rsa, system- 
ic artery resistance, Rsmabp, the resistance of IABP-implanted systemic arteries; Csa, 
systemic artery compliance; Rsv, systemic vein resistance; Csv, systemic vein compli- 
ance; Rtr, tricuspid valve resistance; Cra, right atrium compliance; and Rpu, pulmonary 
valve resistance. Csajabp is calculated as the product of Csa and a scale factor for the 
IABP effects. 



Fig. 2. Electrical activation time mapped to mechanical component of ventricular 
computational mesh. The activation time is defined as the instant at which trans- 
membrane voltage exceeds 0 mV. EAT indicates electrical activation time. 

Table 1. Flow, stroke volume, ejection fraction, pulse efficiency for failure and patho- 
logical model according to the compliance change 





Heart failure considering 


References 


1 


Remodeled geometry 


Gurev et al. (8) 


2 


Reduced electrical conductivity 


Helmet al. (10) 


3 


Increased stiffness of failing myocardium 


Wu et al. (11) 


4 


Systolic dysfunction 
(remodeled cellular Ca transient) 


0'Rourkeet al. (12) 
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the peak and time constant of the calcium transient function 
([Ca](t) of Eq. 55 in the myofilament model of Rice et al. (7)), 
which served as an input, were reduced to 70% of the normal 
values. This resulted in a calcium transient with a reduced peak 
and prolonged relaxation rate, which are the key features of the 
remodeled Ca transient in HF (12). 

Combining the ventricle model with the circulatory and 
IABP models 

We combined the FE heart model with a lumped model of the 
circulatory system and IABP to predict the cardiovascular re- 
sponse during IABP therapy. The model is formulated as an elec- 
trical analog model consisting of elements such as resistors, ca- 



pacitors, and diodes. Using the hemodynamic parameters re- 
lated to blood circulation, this model calculates the pressure, 
volume, and blood flow for each of the six compartments mak- 
ing up the vascular system and atria (Fig. 1): the systemic arteries 
and veins, pulmonary arteries and veins, left atrium, and right 
atrium. Individual elements in every compartment are mod- 
eled by considering the local resistance to blood flow and com- 
pliance of blood vessels. The atria are modeled as two-chamber 
pumps with time-varying compliance. The inertia of blood is 
ignored due to its relatively small effect on the canine cardio- 
vascular response. The four heart valves that surround the ven- 
tricles and control the blood flow direction are modeled by di- 
odes and resistances. 




0 1 1 1 1 1 0 1 1 1 1 1 — 

0 600 1,200 0 600 1,200 

Time(ms) Q Time(ms) Q 



Fig. 3. Simulated pressure waveform in the LV and systemic artery. HF without IABP support (A), and HF with the IABP operating at levels 1 (B), 2 (C), 3 (D), 4 (E), and 5 (F). 



http://dx.doi.org/! 0.3346/jkms.201 3.28.1 .93 



http://jkms.org 95 



JKMS 



Lim KM, et al. • Computational Modeling of Intra-Aortic Balloon Pumping 



Briefly, the IABP component was modeled as the time-vary- 
ing compliance of the systemic arteries. To generalize the pat- 
terns of inflation and deflation of an IABP, a harmonic wave- 
form was used for the time-varying compliance of the systemic 
arteries. The harmonic waveforms for the compliance of the 
aorta were expressed as follows: 



CsA,IABP (t) = CsA X SFlABP (t) 

SFiabp (t) = (1-sf) + sfx cos - 0) 



(1) 

(2) 



Here, Csa,iabp is the time-varying compliance of the systemic 
arteries with the IABP, Csa is the compliance of the systemic ar- 
teries without the IABP, SFiabp is a scale factor for the IABP, sf is 
the level of the scale factor, BCL indicates the basic cycle length 
of the ventricle, and 0 indicates the time shift between the ven- 
tricular contraction cycle and IABP inflation cycle. The sf is ma- 
jor parameter which is proportional to the stroke volume of IABP 
and related to the pumping compliance of IABP. In this study, sf 
was given values of 0.05, 0.1, 0.15, 0.2, and 0.25, and 0 was set to 
3.66 radians (350 ms shift), instead of end-systole. 

RESULTS 

Fig. 2 shows the simulated pressure waveforms in the LV and 
systemic artery under HF (A), and HF with the IABP operating 



at level 1 (sf = 5%, B), 2 (sf = 10%, C), 3 (sf = 15%, D), 4 (sf = 20%, 
E), and 5 (sf = 25%, F). The mean arterial pressure under HF was 
125 mmHg, which matches that in congestive heart failure. As 
the level of IABP treatment was increased, the aortic pressure 
during diastole, which is related to coronary perfusion, increased 
and the left ventricular peak pressure during systole, which is 
related to the after-load, decreased. 

Fig. 3 shows the calculated transmural distribution of the con- 
tractile ATP consumption rates at end systole and end diastole 
for HF without an IABP (Fig. 3A), and the model with the IABP 
operating at level 5 (Fig. 3B). The ATP consumption rate changed 
significantly in the left ventricle rather than the right during sys- 
tole. The contractile ATP consumption decreased significantly 
in the presence of the IABP. Compared with the case involving 
HF without an IABP, the contractile ATP consumption decreased 
by 12% with IABP therapy (see the ATP consumption rate in 
Table 2). 

Fig. 4 shows the calculated transmural distribution of the fi- 
ber strain at end systole and end diastole for HF without an IABP 
(Fig. 4A) and the model with the IABP operating at level 5 (Fig. 
4B). The fiber strain decreased significantly in the presence of 
the IABP. 

Fig. 5 shows the pressure-volume curves for the six cases stud- 
ied: HF without IABP therapy, and HF with the IABP at levels 1 
to 5. The pressure-volume curves shifted to the left and down- 



Table 2. Hemodynamic responses for HF and HF with the IABP according to the compliance change 



Responses 


No IABP 


IABP (level 1) 


IABP (level 2) 


IABP (level 3) 


IABP (level 4) 


IABP (level 5) 


CO (mL/s) 


52.93 


54.11 


55.31 


56.54 


57.8 


59.07 


SV (mL) 


31.76 


32.46 


33.19 


33.93 


34.68 


35.44 


SW (mmHg*L) 


3,495 


3,492 


3,487 


3,481 


3,471 


3,457 


LVPP 


140.6 


137.5 


134.5 


133 


130 


127.5 


MAP 


128 


130 


133 


135 


138 


141 


ATP consumption rate 


319 


311 


303 


296 


288 


281 


EF (%) 


22.75 


23.51 


24.48 


25.48 


26.68 


27.62 


IABP, intra-aortic balloon pump; CO, cardiac output; SV, stroke volume; SW, stroke work; LVPP, left ventricular peak pressure; MAP, mean arterial pressure; EF, ejection fraction. 
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Fig. 4. Transmural distribution of the ATP consumption rate. Heart failure ventricles without IABP support (A) and with IABP support at level 5 (B). 
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Fig. 5. Transmural distribution of the mechanical strain. Heart failure ventricles without IABP support (A) and with IABP support at level 5 (B). 
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Fig. 6. The pressure-volume curves for the six cases studied. Heart failure without 
IABP therapy, and HF with the IABP at levels 1 to 5. 

ward. The end-diastolic volume decreased by 8 mL (Fig. 6); this 
indicates that the IABP results in ventricular volume unloading 
by reducing the end-diastolic volume. 

Table 2 shows the hemodynamic responses for HF and HF 
with the IABP at levels 1 to 5. While the cardiac output and stroke 
volume increased, the stroke work and ATP consumption de- 
creased with IABP therapy due to the pressure unloading effect. 
The decreased left ventricular peak pressure also indicates the 
pressure unloading effect during IABP therapy. 

DISCUSSION 

This study investigated the changes in the contractile energy 
consumption of the myocardium with IABP with different in- 
flation volumes using a cardiac electromechanical model with 
a lumped model of the circulatory system and IABP function. 
To generalize the patterns of inflation and deflation of the IABP, 
a harmonic waveform was used for the time-varying compliance 
of the systemic arteries. According to the IABP inflation volume, 



the arterial pressure increased during diastole, while the ventric- 
ular peak pressure decreased, which indicates improved coro- 
nary perfusion and ventricular unloading. To estimate the change 
in the ventricular workload directly, the contractile ATP con- 
sumption of the ventricle was computed (Fig. 3). With greater 
inflation, the ATP consumption decreased significantly. The 
transmural distribution of the mechanical strain was also com- 
pared between HF ventricles and HF ventricles with an IABP at 
level 5 in order to estimate the volume-unloading effect of the 
IABP (Fig. 4). The fiber strain also decreased significantly with 
an IABP. The pressure-volume diagrams in Fig. 5 were consis- 
tent with the fiber strain distributions. 

The contractile ATP consumption and fiber strain varies trans- 
murally. In general, the sarcomere length affects the density of 
effective cross-bridges and therefore the contractile force; this 
is known as the Frank-Starling law of the heart. A longer sarco- 
mere length results in greater contractile ATP consumption be- 
cause of the greater density of the attached cross bridges in my- 
ocytes, as expressed in Eq. 3 in the Appendix. Previously, we 
demonstrated that myonbers in the endocardium have greater 
end-diastolic lengths than myonbers in the mid-myocardium 
and epicardium (6). Fig. 4 also shows that the end-diastolic strain 
and ATP consumption rate of the endocardial sarcomere are 
greatest. This explains why the contractile ATP consumption by 
the endocardium was greater than that at the mid-wall or epi- 
cardium (Figs. 3 and 4). The cardiac output, stroke volume, stroke 
work, MAP, ATP consumption rate, and ejection fraction are all 
increased with IABP treatment (Table 2). Interesting point is 
that although the ATP consumption is decreased significantly, 
stroke work is decreased only slightly, which indicates that the 
IABP helps the failed ventricle to pump blood efficiently. 

There are several limitations. We used an image-based elec- 
tromechanical model of failing canine ventricles, in which the 
ventricular geometry, fiber architecture, and hemodynamics 
differ from those of human ventricles. We also considered only 
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the contractile ATP consumption of the myocardium, although 
ATP is used in other ways, such as by sarcoplasmic/endoplasmic 
reticulum calcium ATPase (SERCA), plasma membrane Ca- 
ATPase (PMC A), and Na pumps. Their functions are changed 
at failing heart condition (13-15). Finally, we did not considered 
remodeled electrical activities of myocardial cell in failing heart. 
We assumed that the electrical conductivity changes in tissue 
level encapsulated those changes in cellular level. However, 
these limitations should not alter the main findings of this study 
greatly. 

We conclude that LABP improves coronary perfusion as well 
as pressure unloading and volume unloading and also helps 
the failed ventricle to pump blood more efficiently. 
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■ APPENDIX ■ 



Contractile ATP consumption of myocardium 

The contractile energy consumption of the myocardium was quantified by calculating the contractile ATP consump- 
tion in the myofilament model of Rice et al. (7). The contractile ATP consumption rate, E, per unit volume was calculat- 
ed as a function of the ATP- consuming cross-bridge detachment rate {gxbr) and the single overlap fraction of thick fila- 
ments {SOVFmck), 



E = g X bT x SOVFmck 



(3) 



where gxbr, which indicates the ATP-consuming detachment transition rate, and SOVFmck, which indicates the single- 
overlap fraction of the thick filament, were functions taken from the original myofilament model of Rice et al. (7). 
Then, g X br can be derived as: 



TmpC-37 



gxbT = gxb x gxbmd x xbmodspedes x Qg X b ^ 10 



(4) 



The term xbmodspedes scales all cross-bridge cycling rates to account for species-based differences (0.2 for canines). 
The temperature dependence of gxbr (Qgxb) is set to a default Q10 value of 6.25. The term g X br had strain dependence in 
the rate modifier gxbmd, defined as: 



gxbmd - 



(( x 0 - xXBpostR \ 2 \ 
J J 

f (x 0 -xXBpostR\ 2 \ 

6x pH Xo n 



if xXBpostR < Xo 
if xXBpostR > Xo 



(5) 



where a p and a n set the effects of strain for positive and negative shortening velocities, respectively; x Q is the distortion 
of the cross-bridge head; xXB Pos tR is the mean distortion of states XBp 0S tR, which is the probability of the transition to the 
post-rotated force-generating status. The ATP-consuming detachment transition rate g X br is proportional to gxbmd, 
which depends on the strain of the cross bridges, and in turn, this strain depends on the velocity of contraction. See Rice 
et al. (7) for details. Finally, we calculated the ATP consumption distribution by integrating the ATP consumption rate 
with time for 600 ms and the ATP consumption of the entire ventricle by spatially integrating the local ATP consump- 
tion through the entire ventricular volume. 
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